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ABSTRACT 

The rotation curve for the IV galactic quadrant, within the solar circle, is derived from the Columbia 
University - U. de Chile CO(J=1^0) survey of molecular gas. A new sampling, four times denser 
in longitude than in our previous analysis, is used to compute kinematical parameters that require 
derivatives w/r to galactocentric radius; the angular velocity the epicyclic frequency k{R), and 

the parameters A{R) and B{R) describing, respectively, gas shear and vorticity. The face-on surface 
density of molecular gas is computed from the CO data in galactocentric radial bins for the subcentral 
vicinity, the same spectral region used to derive the rotation curve, where the two-fold ambiguity 
in kinematical distances is minimum. The rate of massive star formation per unit area is derived, 
for the same radial bins, from the luminosity of IRAS point-like sources with FIR colors of UC HII 
regions detected in the CS(J=2^1) line. Massive star formation occurs preferentially in three regions 
of high molecular gas density, coincident with lines of sight tangent to spiral arms. The molecular 
gas motion in these arms resembles that of a solid body, characterized by constant angular velocity 
and by low shear and vorticity. The formation of massive stars in the arms follows the Schmidt law, 
r^MSFR (X [Egas]", with an index of n = 1.2 ± 0.2. Our results suggest that the large scale kinematics, 
through shear, regulate global star formation in the Galactic disk. 

Subject headings: Galaxy: kinematics and dynamics — structure, ISM: molecules, stars: formation 


1. INTRODUCTION 

The rotation curve, describing the circular speed of ro¬ 
tating material as a function of galactocentric radius, is 
a fundamental tool for the study of the kinematics of 
our Galaxy. It is best derived, because of interstellar ex¬ 
tinction, from observations of atomic and molecular gas 
in radio and mm wavelengths. The derivation involves 
determining the terminal velocity, or maximum absolute 
radial velocity relative to the Sun, toward lines of sight 
that sample the Galaxy within the solar circle (quad¬ 
rants I and IV). Such terminal velocities correspond, as¬ 
suming pure circular motion, to the tangent points to 
circumferences around the galactic center, named sub¬ 
central points. These points subtend a circumference 
that connects the solar position with the galactic center. 
A detailed analysis of the rotation curve can reveal im¬ 
portant physical characteristics of the rotating material, 
such as the amount of shear and vorticity at each galac¬ 
tocentric radius. These physical quantities regulate the 
gravitational stability of a differentially rotating gaseous 
disk and, consequently, the large scale distribution and 
properties of star formation in the galactic disk. 

The first derivation of the rotation curve for the IV 
galactic quadrant that made use of the GO(J=1^0) line 
- the best tracer of molecular hydrogen in the interstel¬ 
lar medium - was presented by Alvarez, May, & Bronf¬ 
man (1990). The spectral data used to determine the 
terminal velociti es were taken from the Golumbia - U. 
de Ch ile surveys ijGrabelskv et aljfl987l: iBronfman et al.l 
[T^ . which have a sampling interval of 0°.125 (roughly 
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the beam siz e). _ However, the terminal velocities in 

lAlvarez et alJ 11199(111 were measured only every 0°.5 in 
galactic longitude, due to difficulties involved in the vi¬ 
sual examination of a very large number of spectra. A 
new derivation of the rotation curve, that uses a com¬ 
puter search code to examine all the available spectra 
( RilSOOO), is presented here. The disk kinematic char¬ 
acteristics in the IV galactic quadrant are analyzed in 
detail, from this new rotation curve. These characteris¬ 
tics, as a function of galactocentric radius, are compared 
with the molecular gas density and with the local rate of 
massive star formation. 

A proper derivation of the spiral pattern of our Galaxy 
requires knowledge of the distances to the adopted trac¬ 
ers. These distances are also required to compute the 
masses and luminosities of such tracers. For the gas, 
kinematical distances can be obtained from radial veloc¬ 
ity data of radio line observations, adopting a rotation 
curve, under the assumption of pure circular motions. 
For clouds within the solar circle, however, there is a two¬ 
fold ambiguity in the kinematic distance, that is difficult 
to circumvent and has to be resolved in a case-by-case ba¬ 
sis. But in the vicinity of the subcentral points such am¬ 
biguity is minimal, since at the subcentral points them¬ 
selves the kinematic distances are univocally defined. 

It is worth noting that large scale streaming motions in 
spiral arms, with amplitude of ~10 km/s, which produce 
deviations from pure rotation, have been observed in a 
number of regions of the Galaxy (Burton et al. 1988). 
Streaming motions of such amplitude may introduce un¬ 
certainties of up to 5% in the estimation of galactocen¬ 
tric radii when the streaming is along the line of sight. 
In such unfavorable case, the corresponding uncertainties 
in the estimated distances, for the section of the Galaxy 
analyzed here, may go from of 0.6 kpc to 1.7 kpc. In any 
case, for objects beyond ~ 3 kpc from the Sun, because of 
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optical extinction, kinematical distances are usually the 
only ones available. 

Massive stars are formed within aggregates of molec¬ 
ular gas and dust of 10^-10® solar masses, about 50-100 
pc in size, which are commonly known as giant molecu¬ 
lar clouds, or GMCs for short. The association between 
OB stars and the interstellar medium has been estab¬ 
lished through optical, infrared, and CO observations of 
GMCs close enough to be largely unaffected by extinc¬ 
tion (Orion, Carina, etc). The physical conditions in 
GMCs control their rates of OB star formation, and are 
one of the ma in agents tha t regulate the evolution of the 
galactic disk l|Evanslll999|l . 

There is a close relationship between the galactic spi¬ 
ral structure and the formation of GMCs and, hence, 
with the formation rate of OB stars lIDame et al.lll98ft 
LSolomon et alJlT98firi . Therefore, the GMCs and the re¬ 
gions of OB star formation provide a very good tool 
to trace the spiral arm pattern of a galaxy. An early 
descript ion of the Milky Way spiral arm pattern was 
given bv iCeorgelin fc CeorgelinI lll97fill . who observed the 
H109q! line emitted in HII regions associated with young 
massive stars. A four arm spiral pat tern for the south¬ 
ern M ilky Way was later proposed bv iCaswell fc HavnesI 
iTT^ . using a larger observational database of hydro¬ 
gen recombination lines (H109a & HI 10a). The four 
arm spiral pattern is in general agreement with that 
obtained fr om HI and CO large scale observations of 

Star formation is likely to occur in regions where the 
gas in the Galactic disk is unstable to the growt h of grav¬ 
itation al perturbations. In a classical paper, ISchmidll 
(TT^ introduced the parametrization of the volume 
density of star formation and the volume density of gas, 
relating them through a power law; such parametriza¬ 
tion, know n as ’’Schmidt Law”, has been studie d obser- 
vationally llKennicuttlll989(: lWon^^fc^^l^| 2002ll and ex - 
plained on theoretical grounds llToomrell 964 lTa,nll2nnFll . 
A study of the gas stability in the galactic disk must in¬ 
clude (a) comparison of the gas density with a critical 
value above which the gaseous aggregates undergo grav- 
itational collapse llToomrel 11 964t iKenniciitn 11 989fl and 
(b) examination of the gas shear rate, that governs the 
process of destruction of molecular clouds (e.g. Ken¬ 
ney, Carlstrom, & Young 1993; Wong & Blitz 2002), 

J iresumably through the in jection of turbulent motions 
Mac Low fc Klessenll20(^. 

The link between massive star formation and kinemat¬ 
ical conditions in disks has been studied mostly for ex¬ 
terna l spiral galaxie s JM Ito et al.lll999t iWong fc Blitd 
l2002t iBosssier et al.ll2003^ . where the spatial resolution 
that can be achieyed by the observations is not as good 
as for the Milky Way. The main goal of the present paper 
is, therefore, to accurately describe the spiral arm struc¬ 
ture in the subcentral vicinity of our Galaxy, focusing on 
the molecular gas kinematics, density, and on the rate of 
massive star formation, with the hope of contributing to 
the understanding of the formation and evolution of disk 
galaxies in general. The analysis is carried out for the IV 
galac tic quadrant, where th e spiral structure is more evi¬ 
dent l|Bronfma,n et a,1.l198Sjl than in the I quadrant. Pre¬ 
limina ry work has been presentedJjy jCarrascofc Serran^ 
iTT^ and, more recently bv iLuna et al .1 (1299111 . 


Section (§2) describes the obseryational datasets used, 
the most complete ayailable in their kind. These data 
are used in section (§3) to deriye the rotation curve and 
analyze the relation between molecular gas kinematics, 
molecular gas surface density, and massive star formation 
rate. The validity of Schmidt Law for the Milky Way is 
analyzed in section (§4), and a summary of the results is 
given in section §5. 

2. OBSERVATIONS 

The data used to derive the rotation curve and the 
molecular gas surface density are part of the Columbia- 
U. Chile ^^CO(J=1^9) surveys. These surveys pro¬ 
vide us with the most extensive and homogeneous 
observati onal dataset of CO e mission in the galac¬ 
tic disk llGrabel skv et alJ 119871: iBronfman et al.1119891: 
iDa.me et al.ll299Hl~ The beam-size of the antenna in the 
CO line is 8'.8, and an angular sampling of 9°. 125 was 
adopted. The surveys cover the entire IV galactic quad¬ 
rant in longitude, and ±2° in latitude about the galac¬ 
tic equator. The velocity resolution is 1.3kms“^, and 
the typical rms noise antenna temperature of the ob¬ 
servations is 9.1 K. Main beam temperatures, Tmb, are 
used throughout the analysis, obtained by dividing the 
antenna temperature TX by the main beam efficiency, 
Tmb = T'^/riMB, with r]MB = 0.82 (IBronfman et alJ 
Elia- Hereinafter we refer to the ^^CO(J=1^0) line 
as CO. For a detailed d escri ption of the obse r vation s see 
iCrabelskv et al.1 (1198711 and IBronfman et 

The rate of massive star format ion is estimated from 
the integrated FIR luminosity (|Boulanger fc Perauld 
Il988^ of IRA S point-like sources, with F IR colors of UC 
HII regions (I Wood fc Churchwelll Il989ll . detected in a 
CS(J=2—>1) survey of 1427 sources in the whole Galaxy 
by Bronfman, Nyman, & May (1996). The CS(J=2—f-l) 
emission line requires high molecular gas densities, of 10^ 
- 19®cm“^, to become excited. Therefore, it constitutes 
a good tracer of massive star forming regions. The survey 
used here is the most complete currently available, list¬ 
ing 843 massive star forming regions in the galactic disk; 
the observed velocity profiles provide a good estimator 
of kinematical distances to the sources. 

The observations of the CS(J=2^1) line toward the 
IV galactic quadrant were obtained with the SEST 
(Swedish-ESO Sub-millimeter Telescope) at La Silla Ob¬ 
servatory, in Chile, with a beam-size of 59”. Typical 
rms noise in the spectra is 9.1 K, at a velocity reso¬ 
lution of 9.52 kms“^. The radial velocity coverage is 
269kms“^, large enough to detect all the sources in 
the radial velocity range allowed within the solar cir¬ 
cle. Hereinafter we refer to the CS(J=2^1) transition 
as CS , and to the IRAS point-like sources detected 
in CS as IRAS/CS sources. A number of 66 sources 
from the IBronfman et al.l (1199611 survey are used in the 
present analysis. This number was complemented with 
13 sources, undetected in the original CS survey, yet de¬ 
tected in a new CS survey with three times better sen¬ 
sitivity, completed with the SEST telescope, to be pub¬ 
lished elsewhere. 

The present analysis adopts the lAU recommended 
values for the galactocentric distance of the Sun and 
for the solar circular veloci ty, i.e. Rn = 8.5 kpc an d 
Vo = 229kms“^ respectively dkerr fc Lvnden-Belll986|l . 
The analysis excludes galactic longitudes from 359° to 
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Fig. 1.— Terminal velocities in the IV galactic quadrant. Top: 
Measurements obtained every (0.°125) in Galactic longitude. Bot¬ 
tom: galactic latitude where each measurement was obtained. 

360°, because the method used to derive the kinematical 
distance is highly uncertain for that region. Further¬ 
more, the kinematics of the central region of the Milky 
Wav are more comp lex than those of the galactic disk 
(iSawada et alJf^Olll . 

3. ANALYSIS 

3.1. Derivation of the rotation curve 

The rotation curve is derived from terminal velocities 
of GO spectra at each sampled galactic longitude 1. Since 
the emission in the IV quadrant is blue-shifted, the most 
negative velocity is recorded, as well as the galactic lat¬ 
itude b where it occurs (Fi g. 1). To defi ne the terminal 
velocity, following w ork bv ISinhal l|1978ll for HI and by 
lAlvarez et all (ll990D for CO, the half intensity point of 
the blue-shifted side of the emission peak with the most 
negative velocity is selected (Fig. 2). To be considered, 
an emission peak is required to be larger than 4 times 
the noise temperature T^ms- The emission lines thus se¬ 
lected have an average peak temperature of 1.8 ± 0.7 K, 
and an average HWHM of 3 ± 0.8kms“^. 

The rotation curve is obtained from the set of termi¬ 
nal velocities by assigning to every longitude a subcentral 
point of galactocentric distance R = Rq \ sin(Z) | and he¬ 
liocentric distance D = Rq cos{l). The rotational velocity 
is then derived for each subcentral point, and assigned 
to the galactocentric radius R = Rq \ sin(Z) |, under the 
assumption of (a) pure circular motions and (b) differ¬ 
ential rotation with a ngular velocity not g rowing with 
galactocentric radius llSofue fc R,ubinll2001fl . The rota¬ 
tion curve obtai ned, shown in Figure 3 (top), is similar to 
that derived bv lAlvarez et 10(11993), but more detailed; 
the agreement for the longitude range covered in both 
studies is better than ±2kms“^ (rms). A slight differ- 



Fig. 2. — The detail of a typical spectrum (/ = 320.000, 
b = —0.250) showing our selection criterion in velocity and tem¬ 
perature. The portion of the spectral velocity zone analyzed is 
marked with crosses at the measured temperature. For this fea¬ 
ture, the maximum temperature is T = 1.83i^ at u = Q7.Skms~^. 
The terminal velocity value Vterm = —72.2km s~^ is selected at 
the half-maximum temperature in the blue-shifted side of the line, 
and is indicated by the vertical dotted line. A 4cr noise level is 
indicated by the horizontal dashed line. 

ence is apparent in the radial range R/Rq = [0.56,0.62], 
where their analysis yields a systematic shift in velocity 
of 2kms“^. 

The distance Z = i?tan(6) of the subcentral point to 
the galactic plane (Fig. 3 bottom) is also determined, 
for each longitude, providing a simple tool to examine 
the Z distribution of the CO emission as a function of 
galactocentric radius. The re s ults o btained compare well 
with those bv lAlvarez et al ] (11993), which adopted the 
same method, and with those obtained from an axisym- 
metric analysis of the CO emission in the IV quadrant 
llBronfman et al.lll988D . 

3.2. Disk kinematics 

The rotation curve contains vast amount of in- 
formation about the kine matics of the galactic disk 
llBinnev fc Tremainel 1198711 . Among the principal pa¬ 
rameters characterizing the kinematics, derived here, are 
A{R) and B{R), which describe the radial trends of shear 
and vorticity, respectively; these parameters, when eval¬ 
uated at R = Rq, correspond to the well known Oort’s 
constants. 

As observed from a non-inertial coordinate frame, e.g. 
that centered on the Sun (Local Standard of Rest, LSR), 
the balance between Coriolis, centrifugal, and gravita¬ 
tional forces induce non-circular orbits about the galac¬ 
tic center. These motions can be described as small pe¬ 
riodic oscillations superimposed onto circular orbits. In 
a coordinate frame with its origin in the galactic cen¬ 
ter and corotating with the solar orbit, these periodic 
oscillations trace small ellipses called epicycles, with an 
associated epicyclic frequency, k, which roughly accounts 
for deviations from circular motion. The value of k{R) 
depends, at each galactocentric radius, on the angular 
velocity D.{R): 
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Fig. 3. — Rotation curve for the galactic quadrant IV. Top: In 
our present analysis, shown by the thin_Jinej_tl^ longhude res¬ 
olution is four times better than in lAlvarez et al.l i 19901) . shown 
by small squares. The rms difference for positions analyzed in 
both studies is 2 kms“^. Bottom: Our results for the position 
Z of the subcentral point w/r to the galact ic equator (thin line) , 
are fairly consistent with those obtained bv lAlvarez et all II199GI) . 
shown b y small squares, and w ith the emission centroid Zq ob¬ 
tained bv IBronfman et al.l 11988^ through an axisymmetric analysis 
of the full quadrant IV CO dataset, (dashed line). 
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or, using t he Oort’s A and B parame ters in their general 
definition ij'Riuuev fc Tremaiuelll DSll . 



R/Ro 


Fig. 4. — Rotation curve and derived kinematical parameters. 
The vertical lines depict the center of the regions presumably with 
solid-body rotation Top: A box size of 0.5 kpc has been used to 
obtain a boxcar smoothed version (thick line) of the measured ro¬ 
tation curve (thin line). Bottom: The epicycle frequency k(R), the 
angular velocity and the Oort parameters A(R) and B{R), 

derived from the smoothed version of the rotation curve. The la 
errors in A(R) and B(R) are shown with thin lines. The segment 
from 0.83 to 0.97 R/Rq (shown at 12Qkm s~^ in the plot) is the 
radial region used to fit the Oort constants A and B. 


grows with galactocentric radius. The shear and vortic- 
ity present three relative minima, coincident with local 
maxima of the epicyclic frequency k{R)^ at radii 0.39, 
0.47 and 0.73 R/Rq. The parameter A{R), directly pro¬ 
portional to the radial derivative of the angular velocity, 
tends to zero at these radii, so the angular velocity is 
almost constant, a characteristic of solid body rotation. 
The relative maxima in k{R) at these three radii can 
be interpretated also as evidence for solid body rotation; 


= ’ BiR) = -niR)+AiR), ( 2 ) 

it is possible to write the epicyclic frequency as: 

k{R) = ^/-A[B{R)][A{R) - B{R)]. (3) 

A smoothed version of the rotation curve (Fig. 4 top) 
is used in the calculation, in order to provide continu¬ 
ity for the derivative of 0.{R) with respect to R . The 
size of the smoothing box is 0.5 kpc, smaller than the 
typical width of spiral arms. There are three galactocen¬ 
tric radial regions, marked in Fig. 4 top, where V{R) 
increases monotonically, so that the radial derivative is 
positive and roughly constant, i.e. the signature of solid 
body rotation. The kinematical parameters A(i?), B{R), 
k{R), and as functions of galactocentric radius, are 

shown in Figure 4 (bottom). 

The angular velocity, fl{R), the epicyclic frequency, 
k{R), and the parameter A{R), describing the shear, de¬ 
crease on the average with galactocentric radius. On the 
contrary, the parameter B(R), describing the vorticity, 


at a given galactocentric radius, k would be ^/2 times 
higher for solid body rotation than for a flat rotation 
curve V{R) — Vq = constant. 

The results of our analysis are consistent with other 
work in that the parameters A{R) and B(R) at i? = Rq, 
are very close to the values of the Oort’s constants An 
and B p recommended by the lAU llKerr fc Lvnden-Bell 
[Mi). The values derived here are obtained from lin¬ 
ear fits to A(R) and B(R), within the range 0.83 to 
0.97R/Ro, extrapolated to R/Rq = 1. Our results yield 
Aq = 14.9 ± 4 kms“^kpc“^, and Bq = —12.3 ± 4 
kms“^kpc“^, in good agreement with those recom¬ 
mended by the lAU, Aq = 14.5 ± 2 kms“^kpc“^ and 
Bq = —12.0 ± 3 kms“^kpc“^. In a similar man¬ 
ner a value of kq = 35 ± 7kms“^kpc“^ is obtained 
for the epicyclic frequency at i? = Rq, in agreement 
with the value of K n = 36 ± 10 kms“^ kpc“^ given by 
iBinnev fc Tremaiii^ lll987ll . The angular velocity ob¬ 
tained for R = Rq, Oq = 27.2 ± 1 kms“^ kpc“^, is com¬ 
parable with the value of Op = 26 ± 4 kms“^ kpc“^, de¬ 
rived for the spiral pattern angular velocity (Fernandez, 
Figueras, & Torra 2001). Our results therefore suggest 
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that the Sun location is close to the corotation radius 
(see also discussion in Amaral & Lepine 1997). 

3.3. The subcentral vicinity 

To study the possible relationships between the kine- 
matical characteristics of the Galaxy, the molecular gas 
density distribution and the massive star formation rates 
(MSFR), it is necessary to use consistent data sets. The 
densities and MSFRs should be estimated, preferentially, 
for the very same regions for which the kinematics have 
been inferred. The rotation curve is obtained using in¬ 
formation from the subcentral points of the disk, so the 
densities and MSFRs we wish to examine are evaluated 
in regions adjacent to the subcentral points. Such pro¬ 
cedure advantageously minimizes the two-fold ambiguity 
in kinematical distance occurring within the solar circle 
since, at the subce ntral points, kinematical di stances are 
univocally defined llBinnev fc Tremainelll98'^ . 

Limiting the analysis to a reduced area of the 
longitude-velocity space avoids also the azimuthal aver¬ 
aging of physical conditions which pertain both to spiral 
arms and to inter-arm regions, a difficulty that appears 
in axisymmetric analysis of the entire longitude-velocity 
space. Azimuthal averaging, over disk annuli, of the gas 
surface density and star formation rates, allows deriva¬ 
tion oftheatajformationa£^it 2 _thresholdjndisk_galax- 
ies JMartin fc Kennicuttll200lt IBosssier et aUlToO.'ili : the 
detailed information on the spiral structure of the disk, 
however, is washed out in the azimuthal average. 

The subcentral viciniti /’, used to inspect the gas phys¬ 
ical conditions and MSFR, is defined here, at each galac¬ 
tic longitude, as the velocity span between the terminal 
velocity and the terminal velocity plus AR = 15kms“^ 
(Fig. 5). The value of 15kms“^ is large enough to 
include all the emission from each CO profile compo¬ 
nent used to select a terminal velocity. Were the emis¬ 
sion to be interpreted as originated by a set of clouds 
at different velocities, the range Ai? implied by AV = 
15kms“^ would be 0.5 kpc at I = 307° and 0.2 kpc at 
I = 342°, roughly the longitude limits of the analyzed 
region. These values of AR are small compared to the 
large scale trends with respect to galactocentric radius 
and therefore do not affect our conclusions. 

The properties of the ISM, molecular gas face-on sur¬ 
face density and MSFR per unit area, are averaged 
within the subcentral vicinity, to evaluate their depen- 
dance on galactocentric radius. It is assumed that the 
gas shares the kinematics of the corresponding subcen¬ 
tral point at each longitude. The molecular gas face-on 
surface density is computed from the same CO database 
used to derive the rotation curve, while the MSFR is com¬ 
puted from the FIR luminosity of the IRAS/CS sources. 

3.4. Molecular gas surface density 

The derivation of molecular clouds mass from our CO 
survey data is based on the widely made assumption 
that the velocity-integrated CO intensity is directly pro¬ 
portional to the total H 2 column density in molecular 
clouds, and is therefore proportional (with a correction 
to account for hel ium) to the total molecula r gas mass 
column density. l|Binnev fc Tremain3ll987|l . To eval¬ 
uate the molecular gas surface density, Sgasj the CO 
emission is binned in galactocentric annuli of thickness 
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Fig. 5. — The subcentral vicinity borders, separated by AF = 
15/cm are shown with the thick lines atop the longitude-velocity 
diagram of CO emission (gray scale), integrated in latitude over 
—2° < 5 < 2°. Horizontal dotted lines mark the regions with solid 
body rotation, as in Figure 4. 

0.05i?/i?o and, for each radial bin, integrated over all ve¬ 
locities within the subcentral vicinity and over [—2°,2°] 
in latitude. The resultant quantity is divided by the sam¬ 
pled face-on area. The proportionality factor X adopted, 
where N(H 2 ) = AkF(^^CO). is equ al to 1.56 x 10^° cm“^ 
(Kkms“^)“^, liHunter et al.lll99% . A correction of 1.36 
is used to account for helium (see Murphy & May 1991). 
The results obtained for the gas surface density in the 
subcentral vicinity, within the solar circle, are shown in 
Figure 6 {toy). 


3.5. Rate of massive star formation 

Massive stars form in the dense cores of giant molecu¬ 
lar clouds. The UV radiation from the young stars heats 
the surrounding dust, which re-radiates the energy prin¬ 
cipally in the far infrared. Most of these cores are de¬ 
tected as IRAS point-like sources, with FIR colors typ¬ 
ical o f ultra-compact HIT regions l|Wood fc Churchwelll 
IT^ . The density is high enou gh for the excitation o f 
the CS lin e, above 10 "^ -10® cm“® llBronfman et al .11 99611 . 
Following lKenniciitil lll998aH . the rate of massive star for¬ 
mation as a function of galactocentric radius is computed 
by adding the FIR luminosity of the IRAS/CS sources 
that fall within each radial bin. 


MSFRfir = 6.5 X 10-1° 



[MqVt 1 ]. 


(4) 
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Fig. 6.— Molecular gas surface density and integrated FIR 
luminosity from massive star forming regions. Top: The face-on 
surface density for the subcentral vicinity. Bottom: Integrated 
FIR luminosity per unit face-on area for the subcentral vicinity 
com pared with an axisymm etric analysis of the southern galactic 
disk IBronfman et al.l2nnnib The analysis of the subcentral vicinity 
shows structure related with the spiral pattern. 


The FIR luminosity is evalu ated as Lf/a = 
llBoulanger fc Peraultlll98j^ . where D is the subcentral 
distance, 


4 

Ffir = ^ (5) 

i=l 

and Fv(i) are the FIR fluxes of the IRAS/CS point- 
sources in the four IRAS bands, as published in the Point 
Source Catalog. The galactic FIR face-on surface lumi¬ 
nosity of IRAS/CS sources as a function of galactocentric 
radius, for the subcentral vicinity, is presented in Figure 
6 (bottom), compared with the axisym metric analysis of 
the co mplete southern galactic disk by IBronfman et HI 

l200nll . 

The FIR surface luminosity estimated here for massiye 
star forming regions constitutes only a lower limit, since 
the emitting regions powered by embedded massive stars 
can be more extended than the IRAS resolution. But 
deconvolution of the extended FIR continuum emission 
can be much more difficult because of superposition of 
sources along the line of sight. For point-like sources 
there is no such confusion since all the CS prohles mea¬ 
sured have only one velocity component (IBronfman et all 
and, therefore, the derivation of kinematical dis¬ 
tances is straightforward. 

The derived MSFR per unit area (Figure 7c) has rel¬ 
ative maxima at the same galactocentric radial regions 
than the molecular gas surface density (Fig. 7 b). These 
radial regions are characterized by solid-body like rota¬ 
tion, as shown in the smoothed version of the rotation 
curve displayed in (Fig. 7 a). 

3.6. Gravitational disk stability 

The Toomre criterion for disk stability llToomrel 
II964lBinnev fc TremainellI987r . used here, is described 


through the Q stability parameter for gas, defined by 


Q{R) 


acK(R) 

irGEgasiR)' 


( 6 ) 


where c is the velocity dispersion, G is the gravitational 
constant, is the gas surface density and a is a di¬ 
mensionless parameter that accounts for deviations of 
real disks from the idealized Toomre thin disk, single 
fluid model (Tan 2000). When Q < I a gaseous disk 
is gravitationally unstable. The particular case Q = 1 
defines a critical surface density value, governed by the 
Coriolis force represented by k(R); 


t{R) = 


acK(R) 
ttG ' 


(7) 


The Toomre parameter Q can be expressed, hence, as 
Q = Yicrit/F^gas- The constant a is usually estimated 
defining Q = \ wh ere massive star fo r mation ceases along 
disks of galaxies llKenniciittlll998ll iHunter et alJII99^ 
iPisano et al.ll200(T|l . Here a value of a = 0.08 is found by 
defining Q = lin the galactocentric radial range R/Ro = 
0.4 to 0.45, wh ere the MSFR pe r unit area is close to 
zero. Following llKennicu til l1998bll , a value of 6 kms ^ is 
adopted for the velocity dispersion. The dependance of 
Q on galactocentric radius is shown in Figure 7d. 

To analyze the case of cloud destruction governed by 
the shear rate instead of th e Coriolis force. T oomre’s cri¬ 
terion has been modified bv iFImegreenI ()l 99.3^ . The shear 
rate is described by the Oort parameter A(R)] a new pa¬ 
rameter, Qa, that evaluates the survival of a cloud, is in 
this case 


Qa{R) 


aA2.5cA{R) 
TrG^gas(R) ’ 


( 8 ) 


Clouds are sheared to destruction, in differentially rotat¬ 
ing disks, for gas surface densities larger than the critical 
surface density 


S 


A 

crit 


(R) 


2.5aAcA(R) 

5/g 


(9) 


so that Qa = 'Fcrit/^gas- This modified stability cri¬ 
terion takes into account the important role of shear in 
the destruction of GMCs, and can be used to describe 
how shear controls the rate of star formation in galactic 
disks in the regions where the Coriolis force does not. 
The parameter aA is found, as a, to be otA = 0.2. The 
dependance on galactocentric radius of Qa is shown in 
Figure 7e. 

A complementary analysis of t he stability of t he ga s 
clouds to tidal shear is proposed bv IKennev et alJ (1199,311 . 
who defines a critical tidal surface density Tiude- 


^ a,(R)[3A(A-B)]V2 

i^tide(R) - - ^ -, ( 10 ) 

where (Tz is the velocity dispersion in the z-direction of 
the disk. If the gas density is less than Hude, tidal shear 
will rip apart the clouds. The relative importance of 
tidal shear and gravitational stability can be expressed 
by the ratio of ^ude to the critical surface density for 
gravitational instabilities as defined by Toomre for an 
ideal disk. 
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Fig. 7.— Summary of the results from our analysis of the sub¬ 
central vicinity. The regions with solid-body like kinematics are 
shown by vertical dash-dotted lines. Correlations are discussed in 
the text, (a) Smoothed rotation curve. (6) Molecular gas face- 
on surface density (c) Massive star formation rate per unit area 
^MSFR- (d) Gravitational instability parameter Ticrit/'^gas^ (e) 
shear instability parameter {/) Parameter ^tidef^grav 

comparing the relative influences of gravity and tide in the disk sta¬ 
bility. 


s 


grav 


CK 


( 11 ) 


Assuming that c « , the ratio depends almost entirely 

on the shape of the rotation curve, 


-‘tide 


-‘grav 



( 12 ) 


and tells whether tidal or gravitational forces determine 
the disk stability. The dependance of ’Sude/^grav in 
galactocentric radius is shown in Figure 7/. 


bins of extent 0.05i?/i?Oi with the first bin centered at 
R/Ro = 0.275. Panel c shows the MSFR, in units of so¬ 
lar masses per year per square pc (Smsfa)) averaged in 
the same galactocentric radial bins as for b. Uncertain¬ 
ties are estimated from Poisson statistics. Panel d dis¬ 
plays the gravitational stability parameter Q defined by 
Toom re; panel e the Qa parameter defined bv lElmegreenI 
iTT^ for the case when shear dominates over Coriolis 
force; and panel/, finally, the ratio of the tidal to the 
gravitation surface instabilities (eq. [11]). The ratio / is 
estimated from the parameters A{R) and B{R) derived 
from observed rotation curve. 

The gas surface density and the MSFR are higher than 
average, as apparent from Figure 7, in regions where the 
rotation is typical of a solid body, roughly indicated by 
the vertical lines at radii 0.39, 0.47 and 0.73i?/i?o (see 
also Fig. 4). These three regions are very close to the ac¬ 
cepted values for the tangent points of spiral arms iden¬ 
tified in the literature (e.g. Vallee 2002 and references 
therein); the 3kpc arm approximately at 0.36i?/i?Oi the 
Norma spiral arm at 0.51i?/i?o and the Crux spiral arm 
at 0.77R/Ro. Figure 7d tells us that although gravita¬ 
tional instabilities are present in the Galactic disk, the 
disk is self-r egulated in the sense that Q is of order 1 
everywhere However, as shown in Figure 

7e, where the Coriolis force is replaced by shear as the 
principal agent of cloud destruction, the three regions 
coincident with spiral arm tangents are characterized by 
relative minima. This result is reinforced by the three 
relative minima in Figure 7/ where the tidal shear is 
clearly lower than that needed to disrupt the clouds and, 
therefore, allows them to survive enough so as to collapse 
gravitationally and ultimately form massive stars. 

To estimate the influence of the derived quantities and 
parameters in the process of massive star formation, 
we examine the correlation between Tjmsfr (Panel c) 

and /'^gasi '^crit/^QO'Si ^tide j'Bgrav •, ^ and 

K. A Spearman rank-order correlation coefficient test 
llPress et al lima) is adopted; the ranking and signifi¬ 
cance of each correlation are listed in Table 2. Massive 
star formation is clearly correlated with the gas surface 
density, with a confidence level larger than 90%, and 
less clearly correlated with the angular velocity and the 
epicyclic frequency, with a confidence level larger than 
75% and lower than 85%. For the variables /^gas 
and 'Etide/'^gravi which take into account shear, there 
is no correlation with respect to the NULL hypothesis 
that the dependent variables are drawn from a random 
distribution. These results are discussed in the following 
section, searching for a possible scenario to explain the 
observed correlations. 


3.7. Disk stability, gas density, and massive star 
formation 

Some of our most relevant results are displayed in Fig¬ 
ure 7, which we now describe in more detail. The val¬ 
ues plotted and their errors are also listed in Table 1. 
Panel a, depicts the smoothed rotation curve used to de¬ 
rive the Oort parameters A(i?) and B{R). Uncertainties 
in the rotation curve, of ±2kms“^, are dominated by 
the deviation of the measured subcentral velocities from 
the smoothed curve. Panel b presents the molecular gas 
face-on surface density evaluated in galactocentric radial 


4. DISCUSSION 

4.1. Spiral arms in the Southern Galaxy 

There are regions in the rotation curve (Fig. 7a) that 
present solid-body like kinematics; these regions are co¬ 
incident with the loci of known spiral arm tangents. The 
molecular gas surface density is at least a factor of 2 
larger for the arm regions than for the inter-arm regions 
(Fig. 7b). In those regions where the rotation is typical 
of a solid body, the rate of massive star formation per 
unit area, ^msfr, presents relative maxima (Fig. 7c). 
The relative maxima in the molecular gas surface den- 
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sity and MSFR are coincident with relative minima in 
the kinematical parameters and ’Eude/^grav 

Such results, as shown in Figure 6b, cannot be obtained 
alone from an axisymmetric analysis of the Galactic disk. 

In those radial regions where the kinematics resembles 
that of a solid body, i.e., spiral arm tangents, marked 
with vertical lines in Figure 7, V/R = R:!constant and 
—AfB 1. Stability analysis suggests that molecular 
gas in regions characterized by solid b ody rotati on is sta¬ 
ble {'Sgas ~ ^crit) Or self-regulated (lTa,nll2nnnll : GMCs 
are not destroyed by tidal disruption (Ttude ^ ^gas), 
and there is an enhancement in 'Smsfr- In regions 
that deviate from solid body rotation, the star forma¬ 
tion decreases, implying that shearing motions may dis¬ 
rupt/destroy clouds and, therefore, play an important 
role in regulating large scale star formation in the disk. 

The destruction of clouds by shearing motions, there¬ 
fore, is an important parameter controlling the large 
scale star formation in the disk, particularly in the inter¬ 
arm regions where shear dominates gravitational stabil¬ 
ity. High tidal disruption appears to be the agent for 
GMC destruction in the inter-arm regions, inhibiting the 
formation of massive stars, while low tidal disruption fa¬ 
cilitates the accumulation of gas in the arms, increasing 
the gas surface density locally and promoting the higher 
observed MSFR. Regions of solid body rotation should 
present less collisions due to the low shearing motions 
and, therefore, the MSFR could presumably be limited 
by a process different from shear. 

The relation between solid body motion and enhanced 
star formation has been also found for external galax¬ 
ies. A nearby example of a galaxy in which shear in¬ 
stabi l ity may regulate s tar formation is M33 lIGorbellil 
I2?in3l IHever et al.ll2004|) . Using a shear stability crite¬ 
rion, they find that the predicted outer threshold radius 
for star formation is consistent with the observed drop 
in the Ha surface brightness. Another example is the 
central region of N GC3504. The molecu lar gas kinemat¬ 
ics are derived bv iKennev et alJ lll993ll . and compared 
with the rate of star formation obtained from the Ha 
recombination line. The central region of NGC3504 is 
characterized by solid body kinematics and high rate of 
star formation; they propose that the behavior of star 
formation is strongly influenced by the strength of tidal 
shear, which can help control the star form ation rate via 
cloud destruction. For irregular galaxies, iHunter et al.l 
iTT^ notice that in slowly raising rotation curves, char¬ 
acteristic o f this type o f galax ies, the stability parameter 
derived bv lElmegreenI lll993D (eq. 8), provides a good 
criterion to describe the boundary between cloud sur¬ 
vival and disruption. The conclusions of these authors 
are similar to our results for the Southern Milky Way 
arms, which show solid body kinematics, i.e. that the 
SF is strongly influenced by tidal shear, and in partic¬ 
ular that tidal shear controls the rate of star formation 
through cloud destruction. 

Theoretical work in agreement with the present results 
has been presented bv IRoberts & StewartI II1987D . They 
use a galactic disk model where the ISM is simulated by a 
system of particles, representing clouds, which orbit in a 
one arm spiral-perturbed gravitational field, and include 
dissipative cloud-cloud collision. Their conclusion is that 
the distribution of GMGs and star formation is enhanced 
across the full finite width of the spiral arm and is not 


restricted to either the preshock or postshock regions. 

The agreement between their theoretical results and 
our observations implies that the best tracer for molec¬ 
ular gas spiral structure is the kinematical parameter 
that compares gravitational stability and cloud shear de¬ 
struction (Tttide/'^grav, Fig. 7f), and is also the simplest 
one, because it involves only kinematical parameters that 
can be derived directly from the rotation curve. Spiral 
arm regions can be identified, therefore, as those where 
shear disruption of clouds has less influence than grav¬ 
itational instabilities, so that molecular clouds can pile 
up, increasing their mass and evolving, until they col¬ 
lapse to form stars that will thereby increase the mean 
kinetic temperature. Such increment in temperature in 
th e spiral arm region s has been preliminarily measured 
bv ILuna et alJ II2004D . 

4.2. Massive star formation and the Schmidt Law 

There is a good correlation, as shown in section 3.7, 
between the MSFR per unit area (Fig. 7c) and the 
molecular gas face-on surface density (Fig. 7b) in the 
Milky Way. Such kind of correlation, formally known 
as Schmidt law, has been explored mostly for exter¬ 
nal galaxies. In this section we evaluate, for the sub¬ 
central vicinity of the Milky Way disk, four different 
enunciations of t he Schmidt Law, whose simplest form 
iKennicut ti l 19891) is, in our case, 


'^MSFR OC (Egas) 


(13) 


Previous efforts at derivin g the Schmidt Law f or the 
Milky Way include those of iGuibert et in dzH), who 
find values of n between 1.3 and 2.0, based on the distri¬ 
bution of OB ass ociations, H H regio ns. H I, and early 
GO surveys; and llBosssier et al.l2003^ . who using known 
gaseous and stellar radial profiles, find that the Schmidt 
Law is satisfie d for the Milky Way from R = 4 to 15 
kpc. Recently llKrumholz fc McKeell2005j) . in a general 
theory of turbulence-regulated star formation in galaxies, 
predict the star formation rate as a function of galacto- 
centric radius for the Milky Way taking into account the 
fraction of molecular gas in the form of clouds, and enun¬ 
ciate Schmidt Law in the assumption that the clouds are 
virialized and supersonically turbulent. 

A second form of the Schmidt law, which modulates 
the gas density by a ki nematical param eter, in this case 
the angular velocity LI (iKenniciitdlf 9891) . is given by 


^MSFR OC EgasU. (14) 

A third derivation of the Schmidt Law (lTa,nll2nnnl) , can 
be expressed as 


^msfr « S",,U(1 - 0.7/3), (15) 


with 


_ d[ln{Vcirc)] 
d[ln(R)] 


(16) 


where Vdrc is the circular velocity at a particular galac- 
tocentric radius. This version of Schmidt law has the 
advantage that the proportionality constant can be de¬ 
rived from theory, and includes evaluation of the mean 
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free path for cloud-cloud collision, the fraction of colli¬ 
sions that lead to star formation, and the fraction of gas 
transformed into stars in each collision. 

A fourth derivation, of theoretical character, that as¬ 
sumes molecular clouds to be viri alized and supersoni¬ 
cally t urbulent, has been given bv lKriimholz fc McKed 
<1200,^ . and is expressed here as 

^MSFR OC foMc’^gas^^'^^^ (l"^) 

where fcMC is the fraction of gas in the form of molec¬ 
ular clouds 

The four enunciations of Schmidt Law are tested (Fig¬ 
ure 8) using the values of molecular gas face-on surface 
density and of the MSFR per unit area obtained here 
for galactocentric radial bins in the Milky Way. Al¬ 
though most derivations of Schmidt Law include both 
atomic and molecular gas, our calculation s are limited to 
themolecular gas component of the ISM. iWonsf fc Blitd 
ll2f)nl have found, for a set of seven galaxies, that the 
azimuthally averaged star formation rate per unit area 
correlates much better with than with Y^hi- One 
should note that they required a large, uncertain correc¬ 
tion for extinction to derive the rate of star formation us¬ 
ing Ha images; however, such correction does not affect 
their main conclusion, i. e., that considering the total gas 
density, rather than H/ and H 2 separately, may obscure 
underlying physical processes that are essential to star 
formation. The situation may be presumably different 
in galaxies where the ISM is predominantly atomic; how- 
ever, even f o r M33 , with a predominantly atomic ISM, 
iHever et al.l ll2004l) has shown that the star formation 
rate per unit area still correlates better with E^f^ than 
with E///. 

The Schmidt law is roughly satisfied in all four cases 
analyzed, with a power index close to 1. The disper¬ 
sion of the data points is large if all radial bins are taken 
into account; however, the relation becomes much tighter 
if the radial bins that correspond to interarm regions, 
shown in Figure 8 as thin squares, are removed. For the 
simplest relation (Fig. 8a), Egas oc {^msfr)^, the in¬ 
dex value n = 1.2 ± 0.2 obtained is in g ood agreement 
with those derived for extragalactic disks llWong fc Blitd 
| 2003 ). The relation plotted at Figure 8b, taking into 
account the kinematics through the angular velocity H, 
with an i ndex value n = 1.0 ± 0.1, is in very good agre- 
ment with iKennicuttI lll998bf) . Because they derive an in¬ 
dex from disk-averaged values for a set of galaxies, rather 
than the variation of star formation rate with surface 
density within a galaxy, presumably describing quite dif¬ 
ferent processes, our results show an underlaying physical 
link between Eg^g and ^msfr that dominates both local 
and global scales in the process of star formation. The 
relation for the Schmidt law at Figu r e 8 c, dependent on 
the kinematics, as proposed bv iTanI ll200nH . shows a re¬ 
lation with an index value n = 1.4 ± 0.2, consistent with 
theory. An index of 1.2 ±0.1 is obtained when f itting the 
relation derived bv lKrumIu^k_^£Mrd^gd l)200fTfl . using a 
value of fcMC = 0.25 IjDame^^D^^ . 

It is apparent that the simplest form (Fig. 8a) of 
Schmidt law, that normall y describes the a veraged prop¬ 
erties of disks of galaxies llKenniciittlll 989|i . applies also 
fairly well to "Smsfr in a scale range which is much 
smaller, like the radial bins in the Milky Way. The galac- 



Log ,0 (E,„) 


Log,o (E,„n) 



Log,„ [E,. 3 n( 1-0.7/S)] Log,o 


Fig. 8.— The Schmidt law, relating massive star forma¬ 
tion rate and gas density in the subcentral vicinity of the Milky 
Way, obtained with four different methods. Thick squares rep¬ 
resent spiral arm tangent regions and thin squares interarm 
regions. Linear fits to the radial bins corresponding to spi¬ 
ral arm regions yield tight correlations in all four cases (see 
text), (a) Simplest expression: Log{'EMSFR[^Qiy‘^~^P<^~^) vs 
Log{T,gas[MQpc~^). (b) Expression including the angular ve¬ 
locity explicitly, Log('E aasft ( R)\M(7) 'gr~^pc~‘^kms~^kpc~^). (c) 
Expression derived by ITarJ i20QQH . (d) Expression derived by 
IKrumholz Kr. McKed 1200,^1 

tic disk kinematics, through shear, additionally modu¬ 
lates the global star formation in the Galactic disk; de¬ 
stroys the clouds in the interarm regions and permits 
the piling up of gas in the spiral arms, allowing grav¬ 
ity or another agents like compression and turbulence 
(via supernovas) , to collapse the c louds into stars (see 
recent reviews bv lElmegreenl (|2002D ; iMac Low KlessenI 

(l2finl . 


4.3. Timescale for the growth of gravitational 
instabilities 

The gas depletion timescale due to star formation, tsf, 
can be crudely estimated from the derived MSFR per 
unit area, T^msfr, and the molecular gas surface density, 

Tjgas 5 


TSF 


V 

^gas 

T^msfr 


(18) 


The galactocentric trend of tsf is shown in Figure 9 
(top). Its average value is 10^*^ yr, an upper limit since 
Timsfr represents a lower limit. The timescale for 
growth of GMGs through gravitational instabilities can 
be expressed as r ~Q/k l|Larsonl[l988|l . Results here 
yield timescales of the order ^ lO^yr, with relative min¬ 
ima at the location of the spiral a rm ta ngents (Fig. 9 
middle). Following iKennev et a,1.l ill 993jl . the efficiency 
of star formation (e) can be estimated via tsf=£~^t' if 
the timescales for gravitational instabilities (t) are sim¬ 
ilar to the timescale for cloud collapse. The values of 
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Fig. 9. — Comparison between galactocentric radial variations 
of the gas consumption time, the growth of gravitational insta¬ 
bilities, and the star formation efficiency. Vertical dash-dotted 
lines show regions with solid body - like kinematics (see text). 
(Top) Timescale for gas consumption due to star formation. {Mid¬ 
dle) Timescale for gravitational instability growth. {Bottom) Star 
formation efficiency (percentage) using a constant 10® yr life time 
{tob) for OB associations. 

a few tenths of a percent obtained, however, are rather 
low as compared with values of a few percent typically 
reported in the literature. 

The efficiency for massiv e star formation, according to 
iMac Low fc KlessenI 1120041) , can be further expressed as 


e = tob 


^MSFR 

V 

^gas 


tqb 

TSF 


(19) 


Adopting a typical lifetime for an OB star region of 
TOB = 10®yr and our results for T^msfr and Sg^s we ob¬ 
tain values for e of a few percent in the spiral arms, and 
extremely low values for inter-arm regions, (Fig. 9 bot¬ 
tom) . A comparison of the timescales of gas consumption 
and growth of gravitational instabilities, as a function of 
galactocentric radius, suggests that fast star formation is 
present in regions where rapid cloud growth takes place. 


5. CONCLUSIONS 

This work analyzes the correlation between molecular 
gas kinematical properties, molecular gas surface den¬ 
sity, and rate of massive star formation in the IV galactic 
quadrant, using the most complete data bases available. 
The analysis is carried out for galactocentric radial bins 


A. et al. 

0.5 kpc wide, a compromise to avoid arm-interarm con¬ 
fusion while having good statistics. The data used are 
restricted to the subcentral vicinity, to avoid the two-fold 
distance ambiguity within the solar circle. The main con¬ 
clusions are: 

• The rotation curve obtained is similar to that pre¬ 
sented by Alvarez et al. (1990). Since the sampling in 
longitude is 4 times denser, however, it can be used to 
calculate, as a function of galactocentric radius, kinemat¬ 
ical parameters that require radial derivatives. 

• The angular velocity, the epicyclic frequency, 

k{R), and the parameter A{R) describing the gas shear, 
tend to decrease with galactocentric radius; the param¬ 
eter B{R), describing the gas vorticity, tends to grow. 
The values derived for Oort’s constants Aq and Bq, at 
R = i?o, are consistent with those recommended by the 
lAU. 

• Shear and vorticity have relative minima, and the 
epicyclic frequency relative maxima, at radii 0.39, 0.47 
and 0.73 R/Rq, coincident with the known positions of 
spiral arm tangent regions. Near these radii the kine¬ 
matics are characteristic of solid body rotation: A{R), 
proportional to the radial derivative of the angular ve¬ 
locity, tends to zero, so the angular velocity is roughly 
constant. The relative maxima in epicyclic frequency are 
consistent with such scenario since k is times higher 
for solid body rotation than for a flat rotation curve. 

• Differential rotation and shear are weaker for the 
spiral arm regions than for the interarm regions. The 
relative importance of tidal shear w/r to gravitation in 
the stability of the gas for spiral arm regions, where the 
rotation curve resembles that of a solid body, is half than 
for spiral arms, where the rotation curve is nearly flat. 

• Massive star formation occurs in regions of high 
molecular gas density, roughly coincident with the three 
lines of sight tangent to spiral arms. In these arms 
the formation of massive stars follows the Schmidt law, 
^MSFR cx [Egas]”, with an index of u = 1.2 ±0.2. While 
this law is characteristic of spiral galactic disks, here it 
applies to much smaller spatial scale. A modified ver¬ 
sion of Schmidt law, which modulates the gas density by 
the angular velocity, ^msfr cx Egas D, describes better 
the behavior of the gas at this scale, suggesting that the 
kinematics, through shear, regulate global star formation 
in the Galactic disk. 
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TABLE 1 

Values of our derived and measured parameters in the adopted bins at the subcentral vicinity. Our analysis beyond 

0.875i?/i?o IS uncertain and therefore excluded. 


Gal. radii “ 

circ. vel. 

y c 

^qas 

Smsfr{10 ") '' 




R/Ro 

kms~^ 

Mopc-'' 

Moyr-ipc-'' 





±2 






0.325 

212.4 

3.57 ±0.03 

0.23 ±0.03 

0.99 ±0.11 

1.12 ±0.10 

1.03 ±0.12 

0.375 

209.6 

3.93 ±0.03 

0.89 ±0.15 

1.00 ±0.04 

0.43 ±0.10 

0.56 ±0.11 

0.425 

214.2 

2.33 ±0.03 

0.19 ±0.03 

1.07 ±0.14 

1.45 ±0.11 

1.17 ±0.18 

0.475 

215.7 

3.06 ±0.02 

0.75 ±0.07 

1.01 ±0.05 

0.51 ±0.10 

0.61 ±0.11 

0.525 

223.6 

4.12 ±0.02 

1.21 ±0.34 

0.59 ±0.06 

0.58 ±0.06 

0.93 ±0.11 

0.575 

219.3 

1.94 ±0.02 

0.34 ±0.03 

1.09 ±0.11 

1.13 ±0.11 

0.97 ±0.13 

0.625 

218.0 

1.66 ±0.02 

0.28 ±0.04 

1.14 ±0.12 

1.25 ±0.11 

1.00 ±0.15 

0.675 

214.5 

2.19 ±0.01 

0.09 ±0.01 

0.91 ±0.06 

0.65 ±0.09 

0.75 ±0.13 

0.725 

225.3 

3.92 ±0.01 

0.46 ±0.15 

0.63 ±0.01 

0.09 ±0.06 

0.31 ±0.21 

0.775 

236.7 

2.40 ±0.01 

0.18 ±0.04 

0.99 ±0.02 

0.17 ±0.10 

0.34 ±0.20 

0.825 

234.3 

1.68 ±0.01 

0.18 ±0.04 

0.91 ±0.10 

1.03 ±0.09 

1.03 ±0.22 

0.875 

236.1 

0.84 ±0.01 

0.29 ±0.09 

1.69 ±0.20 

2.03 ±0.17 

1.08 ±0.26 

0.925 

233.1 

0.41 ±0.01 

0.30 ±0.06 

- 

- 

- 

0.975 

233.6 

0.33 ±0.01 

0.03 ±0.01 

- 

- 

- 


“Central bin position 

'“Circular velocity from the smoothed rotation curve at central bin position 
'^Averaged molecular gas surface density by bin 
Averaged rate of massive star formation, per unit area, by bin 

®’^’®Stability parameters evaluated using the interpolated version of averaged surface quantities by bin, and the related kinematic parameter 
from the smoothed rotation curve at each position 


TABLE 2 

The Spearman test to explore possible correlations between T,msfr and Sgas, as well as between T.msfr and 
kinematical and instability parameters presented in Table f. 


correlation with ^msfr 

^gas 

^crit/^gas 




wm 

Rank 

0.53 

-0.02 

-0.24 

-0.26 

QSg 

lilESl 

confidence level 


7% 

55% 

60% 



















